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Abstract 


Experimental  and  theoretical  results  are  presented  for  nonlinear  dispersive 
waveguide  modes  in  an  acoustic  duct.  We  report  observations  of  a  localized 
envelope  which  splits  into  two  disturbances  moving  with  two  different  velocities 
of  propagation.  As  a  consequence,  we  predict  that  if  a  signal  is  amplitude 
modulated  at  the  source,  spatial  beating  between  the  two  disturbances  will 
occur,  and  at  periodic  positions  in  space  the  signal  will  become  frequency 
modulated.  These  results  can  have  applications  in  an  all-optical  AM-FM 
conversion  and  in  high  data  rate  fiber  optic  communications. 
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CHAPTER  1.  INTRODUCTION 


Nonlinearity  can  give  rise  to  behavior  that  is  fundamentally  different  from 
linear  behavior.  Solitons  are  the  most  dramatic  example  where  an  exponentially 
localized  wave  of  constant  shape  is  the  result  of  a  stable  balance  between 
nonlinearity,  which  causes  the  wave  to  shock,  and  dispersion,  which  causes  the 
wave  to  spread.  Among  the  most  well  known  examples  of  solitons  are  the 
envelope  solitons.  These  solitons  are  the  result  of  an  instability  whereby  an 
amplitude  modulation  superimposed  on  a  wavetrain  initially  grows  and  eventually 
self-localizes.  As  we  will  see  below,  this  instability  is  realized  if  the  product  of  the 
nonlinear  coefficient  and  the  curvature  of  the  dispersion  curve  (frequency  vs. 
wavenumber)  is  negative. 

When  the  product  of  the  nonlinear  coefficient  and  the  curvature  of  the 
dispersion  curve  is  positive,  modulations  of  a  nonlinear  wave  are  stable.  In  this 
case,  the  most  striking  feature  resulting  from  the  nonlinear  action  is  that  a 
modulation  of  finite  extent  will  split  into  two  disturbances  moving  with  two 
different  velocities  of  propagation  (Whitham,  1974).  This  behavior  should  be 
contrasted  with  linear  theory  where  such  a  disturbance  will  undergo  distortion 
due  to  dispersion,  but  will  not  split.  To  our  knowledge,  no  experimental  evidence 
of  this  effect  has  been  pursued  to  the  present,  with  the  exception  of  this  work. 

We  can  gain  a  physical  understanding  of  amplitude  modulations  of  nonlinear 
dispersive  waves  by  considering  independently  the  effects  of  dispersion  and 
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nonlinearity  on  a  wave.  Although  the  ideas  below  are  very  general,  we  will 
develop  the  conceptual  background  for  acoustics  in  a  waveguide,  as  this  is  the 
system  that  we  use  for  our  experiment. 

A  harmonic  plane  wave  propagating  in  free  space  or  in  a  waveguide  below 
cutoff  satisfies  the  dispersion  relation  w  =  ck  ,  where  (d  is  the  frequency  of  the 
wave,  k  is  the  wavenumber  of  the  wave,  and  c  is  the  speed  of  sound  in  the 
medium.  That  is  the  frequency  is  proportional  to  the  wavenumber  and  the  phase 
velocity  equals  the  slope  dw/dk  of  the  line  in  the  co-k  plane.  In  this  case,  the 
phase  and  group  velocities  are  independent  of  frequency.  Therefore,  a  group  or 
superposition  of  harmonic  waves  with  different  frequencies  all  travel  together 
with  the  same  speed. 

However,  for  acoustic  propagation  in  a  waveguide  above  cutoff  the  dispersion 
relation  has  curvature  (i.e.,  d^w/dk^  sw^^fk)  0).  Hence  the  phase  velocity  is 
frequency  dependent  and  a  group  of  waves  of  different  frequency  -  and  therefore 
speeds  -  will  change  shape  as  it  propagates.  Since  the  motion  of  the  wave 
group  can  no  longer  be  described  by  a  single  velocity,  a  group  velocity  Vg=do)/dk 
is  defined.  The  group  velocity  gives  the  velocity  of  propagation  of  a  group  of 
waves  extending  over  a  frequency  range  small  compared  to  the  center  frequency 
of  the  group. 

To  clarify  the  concepts  of  the  phase  and  group  velocities,  consider  an 
observer  in  a  frame  moving  with  the  phase  velocity.  In  this  frame  an 
unmodulated  wave  is  described  by  a  stationary  pattern  of  crests  and  troughs 
fixed  in  space  and  the  different  points  on  the  wave  correspond  simply  to  points  of 
different  phase  (Fig.  1.1a).  To  an  observer  in  this  frame,  a  modulated  wave 
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would  also  appear  as  a  pattern  of  crests  and  troughs  fixed  in  space  (Fig.  1,1b) 
However  if  (k)  0,  energy  flows  from  regions  of  high  to  low  amplitude  with 
the  group  velocity  and  the  amplitude  of  the  crests  and  troughs  would  vary  in 
time.  Thus  an  observer  fixed  in  a  frame  moving  with  the  phase  velocity  would 
observe  a  propagating  modulation  with  a  relative  speed  |Vp-Vg|. 


a  ) 


b  ) 


Fig.  1.1  a)  Unmodulated  and  b)  modulated  wave  in  a  frame  moving 
with  the  phase  velocity  Vp. 
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For  large  amplitude  acoustic  waves  nonlinear  effects  become  important.  For 
plane  waves  in  a  nondisperst'.  medium,  distortion  of  the  sound  wave  results 
from  amplitude  dependent  effects.  Since  the  wave  speed  is  proportional  to  the 
square  root  of  the  temperature,  it  follows  that  compressions,  where  the 
temperature  is  higher,  propagate  faster  than  rarefactions,  where  the  temperature 
is  lower.  Also,  at  the  compression  the  particle  velocity  is  in  the  direction  of 
propagation  and  at  a  rarefaction  the  particle  velocity  is  in  the  opposite  direction. 
Shock  formation  of  the  waveform  results  in  the  sawtooth  shape  shown  in  Fig. 

1 .2a  as  the  crests  of  the  wave  overtake  the  troughs. 

For  a  dispersive  medium  shocks  are  not  produced  because  frequency  and 
wavenumber  are  not  proportional  and  thus  new  propagating  modes  are  not 
generated  by  nonlinear  interactions.  However,  due  to  self  interaction  effects,  the 
frequency  of  the  wave  is  amplitude  dependent  and  in  the  weakly  nonlinear 
regime  it  is  of  the  form  w  =  (Oo(k)  +  o)2(k)p^.  In  this  expression,  tOo(k)  is  the  linear 
dispersion  relation,  co2(k)  is  the  nonlinear  coefficient,  and  p  is  acoustic  pressure. 
For  the  case  of  fixed  frequency  and  wj  >  0,  the  effect  of  nonlinearity  is  to 
decrease  the  wavenumber  k  as  the  amplitude  of  the  acoustic  field  increases. 
Equivalently,  for  fixed  frequency,  as  the  amplitude  of  the  wave  is  increased  the 
phase  velocity  increases. 

Because  the  changes  in  the  wavenumber  vary  linearly  with  p^,  the  relative 
phase  between  two  points  fixed  in  space  should  also  vary  in  this  manner.  In 
Chapter  5  we  present  experimental  results  that  verify  this  effect  for  the  (1,0) 
cylindrical  waveguide  mode.  This  result  is  the  first  experimental  observation  of 
the  dependence  of  the  relative  phase  on  the  wave  amplitude. 
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Fig.  1.2 


a) 


before 


after 

b) 


0)  =  COo(k) 


C02(k)p^ 


CO2>0 


a)  Shock  formation  resulting  from  nonlinearity  in  a 
nondispersive  medium  as  viewed  in  a  frame  moving  with 
the  speed  of  sound,  b)  For  co2>0  and  fixed  frequency,  the 
effect  of  nonlinearity  in  dispersive  acoustics  is  to  increase 
the  wavelength. 
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To  understand  the  combined  effects  of  dispersion  and  nonlinearity,  consider 
an  initial  state  of  a  modulated  wave  observed  in  a  frame  moving  with  the  group 
velocity  (Fig.  1 .3a).  Because  of  dispersion,  the  group  and  phase  velocities  are 
different  and  an  observer  in  this  frame  would  see  the  crests  of  the  wave 
propagate.  If  the  phase  velocity  is  bigger  than  the  group  velocity  (as  is  the  case 
for  a  mode  above  cutoff  in  a  waveguide)  the  crests  of  the  wave  will  propagate  in 
the  positive  x  direction. 

Without  loss  of  generality  we  can  assume  that  the  nonlinear  coefficient  co2>0. 
Thus,  the  phase  velocity  in  the  regions  of  higher  amplitude  is  greater  than  in 
regions  of  lower  amplitude.  As  a  result,  an  observer  in  the  frame  moving  with  the 
linear  group  velocity  would  observe  bunching  of  the  crests  at  the  leading  edge  of 
the  modulation  and  anti-bunching  at  the  trailing  edge  (Fig.  1.3b).  For  negative 
dispersion,  co''(k)<0,  the  long  waves  travel  faster  than  the  short  waves  and  an 
instability  sets  in.  Dispersive  effects  cause  the  energy  to  approach  the  crests  of 
the  envelope,  and  the  amplitude  at  the  crest  increases.  This  in  turn  accelerates 
the  instability.  Deep  gravity  waves  on  the  surface  of  a  liquid  is  an  example  of  a 
wave  system  with  negative  dispersion.  For  positive  dispersion,  co"(k)>0,  short 
waves  travel  faster  than  long  waves.  In  this  case  dispersive  effects  cause  the 
energy  to  approach  the  troughs  of  the  envelope  and  the  modulation  propagates. 
In  other  words,  to  an  observer  moving  with  the  linear  group  velocity,  the 
modulation  is  no  longer  stationary.  Acoustic  propagation  in  a  waveguide  above 
the  cutoff  frequency  is  an  example  of  a  wave  system  with  positive  dispersion. 

In  general,  if  the  product  CO2C0"  is  negative  a  modulation  will  initially  grow  and 
the  instability  will  eventually  lead  to  the  formation  of  envelope  solitons.  If  the 
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product  is  positive,  the  modulations  are  stable.  In  Chapter  4  we  calculate 
the  nonlinear  coefficient  for  the  (1,0)  mode  in  a  cylindrical  waveguide  and  show 
that  it  changes  sign  from  positive  to  negative  as  the  frequency  is  increased.  For 
this  mode,  finite  extent  modulations  of  a  uniform  wavetrain  are  stable  below  a 
frequency  approximately  twice  the  cutoff  frequency  and  are  unstable  beyond  this 
frequency . 
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b  ) 


Fig.  1.3  a)  Initial  state  of  a  modulated  wave  in  the  frame  moving 
with  the  group  velocity  Vg.  b)  For  o)2>0,  in  a  frame  moving 
with  Vg  nonlinearity  results  in  bunching  of  the  crests  at  the 
leading  edge  of  the  modulation  and  anti-bunching  of  crests 
at  the  trailing  edge. 
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The  stability  of  the  modulation  has  two  important  consequences.  One  is 
AM-FM  conversion.  The  other  is  pulse  splitting.  Both  effects  can  be  physically 
understood  by  considering  again  a  modulated  wave  in  a  frame  moving  with  the 
linear  group  velocity.  In  this  frame,  an  observer  at  a  fixed  location  would  see 
alternating  bunching  and  anti-bunching  of  the  wave  crests  (Fig.  1.4).  Because 
dispersive  effects  cause  the  energy  to  approach  the  troughs  of  the  envelope,  for 
this  observer  it  would  appear  that  at  some  time  a  frequency  modulation  is 
imposed  upon  the  amplitude  modulated  signal.  The  combined  dispersion  and 
nonlinearity  cause  the  energy  flow  from  the  short  wavelength  to  the  long 
wavelength  of  the  FM  modulation  and  an  amplitude  modulation  180°  out  of 
phase  with  respect  to  the  original  signal  results  at  a  later  time.  In  the  frame 
moving  with  the  group  velocity,  the  process  repeats  periodically  and  an  observer 
in  this  frame  sees  that  the  modulations  experience  beats.  In  the  laboratory 
frame  if  a  source  is  generating  an  amplitude  modulated  signal,  some  distance 
away  it  will  become  frequency  modulated.  This  spatial  beating  of  the  modulation 
can  be  only  be  understood  if  for  each  modulation  frequency  there  correspond 
two  group  velocities  of  propagation. 

The  splitting  of  a  single  high  intensity  modulated  wave  has  been  verified 
experimentally  and  is  reported  in  Chapter  5.  A  single  pulse  modulation  on  a  CW 
is  observed  to  split  into  two  pulses.  The  time  difference  between  them  increases 
both  with  distance  and  amplitude. 
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The  organization  of  this  work  is  as  follows.  In  Chapter  2  we  derive  the 
coupled  wavenumber  and  energy  conservation  equations.  Chapter  3  deals  with 
the  quantitative  derivation  of  group  velocity  splitting.  As  a  consequence  of  group 
velocity  splitting  we  predict  the  existence  of  AM-FM  conversion  in  modulationally 
stable  nonlinear  dispersive  systems.  In  Chapter  4  we  apply  these  results  to  the 
specific  case  of  the  (1 ,0)  mode  in  a  waveguide  of  circular  cross  section.  A 
calculation  of  the  nonlinear  coefficient  coj  for  this  mode  is  presented  and  the 
experimental  consequences  of  these  results  are  analyzed.  In  Chapter  5  we 
present  the  experimental  results  for  a  nonlinear  (1,0)  mode  in  an  acoustic  duct 
for  the  frequency  regime  that  leads  to  modulational  stability.  Finally,  Chapter  6 
summarizes  the  results  of  our  theoretical  and  experimental  investigations  of  this 
phenomena. 
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CHAPTER  2.  ONE  DIMENSIONAL  MODULATION  THEORY 


In  this  chapter  we  review  the  theory  of  one  dimensional  modulation  of  linear 
waves  in  a  dispersive  medium.  We  then  expand  our  discussion  to  include  the 
effects  of  nonlinearity  and  the  combined  effects  of  both  nonlinearity  and 
dispersion.  Our  treatment  parallels  the  work  done  by  Whitham  (1974). 


1.  LINEAR  MODULATION  THEORY  IN  ONE  DIMENSION 
a)  Kinematic  Considerations 

If  the  amplitude  and  propagation  wavenumber  of  a  traveling  wave  vary 
only  slightly  over  distances  of  the  order  of  one  wavelength,  the  wave  can  be 
described  locally  by  assuming  a  slowly  varying  amplitude  and  phase.  A  wave 
field  T  can  then  be  described  by  the  expression 

T  =  a(x,t)e'®‘*'^  ,  (2.1.1) 

where  a(x,t)  is  the  local  value  of  the  amplitude.  Over  small  regions  of  space  and 
short  intervals  of  time  the  phase  0  (x,t)  can  be  expanded  in  a  Taylor  series  and 
up  to  first  order,  we  obtain 

0  =  eo+x— +  t— ,  (2.1.2) 

ox  at 

where  0^  is  the  phase  evaluated  at  the  origin  at  time  t=0. 
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Over  small  regions  of  space  and  during  short  time  intervals,  we  define  the 


local  wave  number  and  frequency  at  each  point  as 


....  50  50 

k(x.t)  =  —  ,  oi(x.t)  =  -— . 

5x  5t 

Thus,  points  of  constant  phase  move  with  velocity 

dx  w 
^  "  'k  ■ 


(2.1.3) 


(2.1.4) 


An  observer  following  a  particular  crest  moves  with  the  local  phase  velocity  but  in 
general  sees  the  local  wave  number  and  frequency  changing:  that  is  neighboring 
crests  get  farther  away. 

By  eliminating  the  phase  0  from  each  expression  in  Eq.  (2.1.3),  we  obtain 

^  +  ^  =  0  .  (2.1.5) 

5t  5t 


Equation  (2.1.5)  has  the  form  of  a  conservation  law,  where  the 
wavenumber  k  in  the  first  term  corresponds  to  the  density  of  waves  or  number  of 
wave  crests  per  unit  length  and  the  frequency  co  in  the  second  term  corresponds 
to  the  wave  flux  or  number  of  wave  crests  crossing  the  position  x  per  unit  time. 

Equations  (2.1.1)  -  (2.1.3)  can  be  extended  to  include  the  case  of  two  and 
three  dimensions.  However,  the  one  dimensional  equations  are  sufficient  to 
describe  waveguide  propagation  and  we  will  use  them  exclusively  throughout  our 
discussion. 
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b)  Energy  Propagation 

For  a  general  dispersion  relation  co  =  co(k)  the  phase  velocity  is  in  general 
frequency  dependent.  Thus,  if  a  disturbance  consists  of  a  group  of  individual 
harmonic  waves  of  different  frequencies  and  hence  speeds,  the  resulting  wave 
shape  will  change  with  position  and  the  motion  of  the  wave  group  cannot  be 
described  in  terms  of  a  single  wave  velocity. 

The  general  form  for  a  propagating  wave  in  a  linear  dispersive  medium 
with  dispersion  relation  co(k)  is 

'F(x,t)  =  |F(k)e''‘’‘-*dk  .  (2.1.6) 


For  a  modulated  wavetrain  with  most  of  the  energy  in  wavenumbers  close  to 
some  value  k^  F(k)  is  concentrated  near  k=ko  and  (2.1.6)  may  be  approximated 
by 


H'(x,t)  =  jF(k)exp[(ikx-i{(o(k,)  +  (k-k,)B'(k„)}l]dk  .  (2.1.7) 

This  in  turn  may  be  written  as  4^  =  a(x-(D'(ko)t)exp{i(koX-(D(ko)t}  where 
a(x-co'(ko)t)  =  jF(ko +K)exp[iK(x-(o'(ko)t]dK  ,  (2.1.8) 


and  where  we  have  substituted  k  =  k^  +  k.  The  function  a(x-(o'(ko)t)  describes 
the  modulations  of  the  wave  field  about  k^.  It  varies  with  x  and  t  in  such  a  way 
that  it  remains  constant  when  viewed  from  a  coordinate  system  moving  with  the 
group  velocity 


do 


(2.1.9) 
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evaluated  at  the  center  frequency  of  the  group,  i.  e.  at  w  =o)(ko).  Thus,  we  say 
that  the  modulation  is  constant  along  the  group  line  x=(o'(ko)  t.  An  observer 
moving  at  the  group  velocity  w'(ko)  always  sees  waves  with  wavenumber  k^  and 
frequency  (o(ko)  but  crests  keep  passing  him. 

The  energy  density  of  a  wave  is  proportional  to  the  square  of  the  wave 
amplitude.  Hence  if  a(x-M'(ko)t)  is  constant  in  a  coordinate  system  moving 
with  the  group  velocity,  the  energy  between  any  two  group  lines  remains 
constant.  In  other  words, 

5  + 1.(036)  =  0  ,  (2.1.10) 

dt  dx  ^ 

where  e  is  the  energy  density.  In  order  to  see  that  (2.1 .10)  is  valid,  consider  the 
quantity 

E(t)  =  pa2dx  =  2Kf^F(k)F(k)*dk  ,  (2.1.11) 

JX|  Jk| 


where  k^  and  kj  are  defined  by 

Xi  =co'(ki)t  ,  X2  =(o'(k2)t  .  (2.1.12) 

If  ki  and  k2  are  held  fixed  as  t  varies,  E(t)  remains  constant.  According  to 
(2.1.12),  the  points  x,  and  X2  are  moving  with  the  corresponding  group  velocities. 
Thus  the  total  amount  of  a^  between  any  pair  of  group  lines  remains  the  same. 
Thus  a^  propagates  with  the  group  velocity.  On  the  other  hand,  e  =  f(k)  a^,  thus 
(2.1.10)  holds  provided 
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The  first  bracket  vanishes  because  E(t)  is  constant  between  any  pair  of  group 
lines.  The  second  bracket  vanishes  in  view  of  (2.1 .5),  the  conservation  of  the 
number  of  waves  crests! 


2.  NONLINEAR  WAVES 

The  modulations  of  a  linear  wavetrain  propagating  in  one  dimension  in  a 
uniform  medium  can  be  described  by  Eqs.  (2.1.5)  &  (2.1.10).  More  generally, 
nonlinearities  modify  the  wave  propagation.  Nonlinearity  results  in  the 
appearance  of  an  amplitude  dependence  the  frequency.  For  small 
amplitudes,  the  frequency  dependence  on  amplitude  can  be  written  as 

(o(k,a^)  =  coo(k)  +  co2(k)a^+...  ,  (2.2,1) 

where  a  is  the  amplitude,  o)o(k)  is  the  linear  dispersion  relation,  and  (02(k)  is  the 
nonlinear  coefficient.  Since  the  energy  e  of  the  wavetrain  is  proportional  to  the 
amplitude  squared,  Eq.  (2.2.1)  can  also  be  written 

co(k,e)  =  too(k)  +  ci)2(k)e+... .  (2.2.2) 

There  are  three  conditions  possible  for  the  nonlinear  coefficient  C02  .  If 
co2(k)  =  0,  the  dispersion  relation  reduces  to  the  linear  case.  If  o)2(k)<  0  both 
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the  frequency  and  phase  velocity  decrease  with  increasing  amplitude.  Finally,  if 
0)2 (k)  >  0,  the  frequency  and  phase  velocity  increase  with  increasing  amplitude. 

An  important  effect  of  the  nonlinearity  is  to  couple  Eqs.  (2.1.5)  and  (2.1.10). 
Thus  substituting  (2.2.2)  into  (2.1.5)  we  obtain 

^  +  (a)i(k)  +  cO2(k)e)^  +  a.2(k)^  =  0  .  (2.2.3) 

For  slow  modulations  of  the  amplitude,  the  last  term  of  Eq.  (2.2.3)  can  be  as 
significant  as  the  term  o)o(k)5k/5x.  However,  the  (o^{k)edk/dx  term  can  be 

ignored  for  weakly  nonlinear  fields  because  it  only  gives  a  small  correction  to  the 
group  velocity. 

Similarly,  in  Eq.  (2.1.10),  we  can  retain  only  the  linear  group  velocity  and 
obtain  the  equations  (Whitham,  1974) 


6k  /  /,  s  ^k  /,  X  ^ 

— +  Qo(k)  — +  o)2(k)— =  0  , 

dt  ox  dx 

(2.2.4) 

de  / .  6e  //  /■  \  5k  . 

— +  coo(k)  — +  a)o  (k)e— -  0  . 
dt  dx  dx 

(2.2.5) 

Equations  (2.2.4)  and  (2.2.5)  constitute  the  basic  description  for  one  dimensional 
nonlinear  waves  in  a  uniform  dispersive  medium.  The  crucial  qualitative  change 
of  nonlinearity  which  couples  both  equations  will  be  explained  in  the  next 
chapter. 
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CHAPTER  3.  GROUP  VELOCITY  SPLITTING 


In  this  chapter,  we  consider  the  general  solutions  to  Eqs.  (2.2.4)  and  (2.2.5). 
Based  on  Whitham's  modulation  theory  we  will  give  a  quantitative  foundation  for 
the  splitting  of  a  modulation  of  finite  extent  into  two  separate  disturbances.  Next, 
in  an  extension  of  these  results,  we  solve  the  eigenvalue  problem  for  a 
monochromatic  space-time  modulation  of  wavenumber  and  amplitude.  From  this 
we  determine  that  the  modulational  energy  is  equipartitioned  Into  two 
propagating  modes.  Lastly,  we  predict  that  a  consequence  of  nonlinearity  on  an 
amplitude  modulated  signal  is  to  introduce  a  frequency  modulation  into  the  signal 
and  that  the  reciprocal  effect  is  also  possible. 


1.  PULSE  SPLITTING 

The  combined  effects  of  dispersion  and  nonlinearity  on  an  initial  disturbance 
or  a  modulation  is  to  couple  Eqs.  (2.2.4)  and  (2.2.5) 

This  system  of  equations  is  hyperbolic  and  in  order  to  determine  the 
characteristic  velocity  for  the  system  (2.2.4)  -  (2.2.5),  i.e.  the  velocity  on  the 
characteristic  along  which  the  wave  propagates,  we  consider  the  linear 
combination* 


*  For  an  excellent  treatment  on  the  classification  of  partial  differential  equations  and  the 
determination  of  the  family  of  characteristic  curves  for  hyperbolic  equations,  the  reader  must 
consult  Whitham's  Linear  and  nonlinear  waves . 
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— +  coo(k)  — +  co2(k)—  +I2  ®  ’  (^•^•^) 

at  ox  ex)  \d\  ax  ax 


which  takes  the  characteristic  form 


dk  dx  dk 


+  U 


( de  dxde 


-0 


(3.1.2) 


provided 


l,^  =  l,Mi(k)  +  l2toi'(l<)e, 
dt 

ri  V 

l2^  =  ljui{k)  +  l,<o;(k). 
dt 


(3.1.3) 


The  system  of  linear  algebraic  equations  (3.1.3)  possesses  a  non-trivial  solution 
when 

^  =  c)[,(k)±(<02(k)coi'(k)e)^"  ,  (3.1.4) 

dt 

and  the  solutions  are  real  (i.e.  propagation  is  possible)  provided  o, (k)^^  (k)  >  0. 

Thus  an  initial  disturbance  or  a  modulation  of  a  source  will  introduce 
disturbances  into  both  families  of  characteristics  and  consequently  two  signals 
moving  with  different  velocities  will  be  detected.  This  behavior  should  be 
contrasted  with  the  linear  theory  where  such  a  disturbance  will  undergo  distortion 
due  to  dispersion,  but  will  not  split. 
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2.  THE  EIGENVALUE  PROBLEM 

In  the  previous  section  we  determined  that  for  co2(k)t0o  (k)  >  0,  two  real 
propagating  solutions  to  Eqs.  (3.1.3)  exist.  In  this  section,  we  approach  the 
solution  of  the  coupled  equations  (2.2.4)  -  (2.2.5)  by  considering  the  condition  on 
the  values  of  possible  frequencies. 

Consider  a  wavetrain  which  has  been  given  a  space-time  modulation  of 
wavenumber  and  amplitude; 

k(x,t)  =  ko +a(x,t)  ,  (3.2.1) 

e(x,t)  =  eo +e(x,t)  ,  (3.2.2) 

where  ko  is  the  unmodulated  wavenumber,  a(x,t)  is  the  modulation  of  the 
wavenumber,  eo  is  the  constant  energy  of  the  unmodulated  carrier,  and  8(x,t)  is 
the  energy  of  the  modulation. 

A  linearization  of  Eqs.  (2.2.4)  and  (2.2.5)  about  the  uniform  state  k^  and  eo 
yields  the  modulation  equations 


da  I  8a  ^8  ^ 

- +  03  n - +  0)  7  —  =  0  , 

at  °  ax  ^  ax 


ae  I  dz 


at 


+  03 


ax 


+  03oeo 


da 

ax 


=  0  . 


(3.2.3) 

(3.2.4) 


If  we  seek  solutions  to  Eqs.  (3.2.3)  -  (3.2.4)  of  the  form  a  =  aoe'™  and 
e  80®'^*"'’*  ■  we  obtain 


21 


(3.2.5) 


(Q-WoKjtto  -C02Keo  =  0  , 

-wdeoKao +(t^-oJoK)eQ  =  0  (3.2.6) 

Equations  (3.2.5)  -  (3.2.6)  have  a  non-trivial  solution  only  if 

Q  =  WqK  ± K  .  (3.2.7) 

The  dispersion  relation  (3.2.7)  shows  that  for  fixed  Q  ,  there  are  two  possible 
wavelengths  of  the  modulated  state  i.e. 


(3.2.8) 


(3.2.9) 


Because  the  nonlinearity  is  weak  (i.e.,  .^codco2eo  «®o)-  (3.2.8)  and  (3.2.9) 
can  be  rewritten  as 


(3.2.10) 


(3.2.11) 
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Because  there  are  two  eigenvalues  to  the  set  of  Eqs.  (3.2.3)  -  (3.2.4),  the 
most  general  solution  for  a(x,t)  and  e(x.t)  for  fixed  must  be  of  the  form 

a(x,t)  =  ,  (3.2.12) 

e(x, t)  =  +  Vje'^-'P*  .  (3.2. 1 3) 


Substituting  solutions  (3.2.12)  -  (3.2.13)  into  (3.2.3)  -  (3.2.4)  and  using  the 
fact  that  the  modes  are  linearly  independent,  we  obtain 

(-Q  +  Wo(k)q)ui  +(02(k)qU2  =  0  (3.2.14) 

(-Q  +  (i)o(k)p)vi  +co2(k)pv2  =  0  .  (3.2.15) 

Using  (3.2.1 0)  and  (3.2.1 1 )  we  obtain 


U2 


(£>2 

Wo  eg 


(3.2.14a) 


V2 


0 


(3.2.15a) 


Equations  (3.2.14a)  -  (3.2.15a)  determine  the  ratios  of  u^  to  Uj  and  to  V2 
respectively.  The  ratio  of  u^  to  v^  and  Uj  to  V2  are  determined  by  boundary 
conditions  on  the  signal. 
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3.  AMPLITUDE  MODULTION  TO  FREQUENCY  MODULATION 

Consider  an  arbitrary  amplitude  modulated  signal  at  x=0  of  the  form 

a(0,t)  =  A[l  +  mcos(Qt)]cos(o)t)  ,  (3.3.1) 


where  A  is  the  amplitude,  m  is  the  modulation  amplitude,  Q  is  the  modulation 
frequency  and  co  is  the  carrier  frequency.  At  the  boundary,  Eq.  (3.2.13) 
becomes 

s(0,t)  =  U2e"'*+V2e'''‘ .  (3.3.2) 

Also,  at  x=0,  the  general  expression  for  the  phase  becomes 

+  (3.3.3) 

q  p 


where  we  have  integrated  the  expression,  — — 

dx 

wavenumber  and  taken  the  associated  arbitrary  function  of  time  in  the  result  to 
be  -cot . 

Because  at  x=0  the  amplitude  modulated  signal  is  of  the  form  (3.3.1),  it 
follows  from  Eq.  (3.3.3) 


t) 


=  ko  +a(x,t)  for  the  modulated 


=0  .  (3.3.4) 

q  P 


This  requires 
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(3.3.5) 


Because  the  nonlinearity  is  weak,  p  is  approximately  equal  to  q  and 

u,  =-v,  ,  (3.3.6) 

which  expresses  the  fact  that  energy  is  equipartitioned  between  the  two 
propagating  modes. 

At  any  point  x,  the  energy  e(x,t)  is 

£(x,t)  =  U2[cos(Qt-qx)  +  cos(Qt-'px)]  ,  (3.3.7) 

where  we  have  considered  the  real  part  only.  Using  trigonometric  identities  and 

defining  A  =  ^  and  ti  =  ,  we  obtain 

2  2 

s(x,t)  =  2u2CosAxcos(Qt-Tix)  .  (3.3.8) 

If  the  modulation  amplitude  m  is  small  compared  to  unity,  from  Eq.  (3.3.7) 
e(0,t)  =  2A2mcosQt  ,  (3.3.9) 

Upon  comparison  with  Eq.  (?  3.2)  we  obtain 

U2=A^m.  (3.3.10) 

Similarly,  at  any  point  x,  the  phase  0(x,t)  is 
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(3.3.11) 


2Ui 

0(x,t)  =  koX-cot - i-sin(Ax)cos(Qt-Tix)  , 

q 


and  from  Eq.  (3.2.14a)  u,  = 


“2  a2 


loi'o'eQ 


A^m 


Also,  using  the  fact  that  q  =  Q  /  aj,  we  obtain  the  following  expression  for  the 
amplitude  modulated  wave; 

a(x,t)  =  A[1  +  2mcos(Ax)cos(Qt-Tix)] 


X  cos 


koX-cot- 


2A^m((0o)^  w 


Q 


'^’060 


sin(Ax)cos(Qt-Tix) 


(3.3.12) 


Thus  a  consequence  of  the  nonlinearity  on  an  amplitude  modulated  signal  is  to 
introduce  frequency  modulation. 

Equation  (3.3.12)  exhibits  reciprocity  regarding  frequency  modulation  to 
amplitude  modulation.  When  cos  Ax  =  0  the  signal  is  entirely  frequency 
modulated.  If  we  consider  this  point  as  the  source  of  the  modulated  wave,  it  must 
follow  that  a  consequence  of  nonlinearity  on  an  frequency  modulated  signal  is  to 
introduce  an  amplitude  modulation  into  the  signal. 


26 


CHAPTER  4.  SOUND  WAVES  IN  A  CYLINDRICAL  WAVEGUIDE 


In  this  chapter  we  review  the  theory  of  sound  wave  propagation  in  a 
cylindrical  waveguide.  Following  Morse  and  Ingard  (1968),  we  examine  the 
general  linear  theory  and  derive  the  waveguide  dispersion  relation.  We  discuss 
the  perturbation  technique  of  Keller  and  Millman  (1972)  that  can  be  used  to 
determine  the  effect  of  small  finite  amplitude  of  an  acoustic  wave.  We  specialize 
the  results  of  the  perturbation  method  to  determine  the  nonlinear  coefficient  toj 
for  the  specific  case  of  the  first  non-axialiy  symmetric  mode  in  a  cylindrical 
waveguide.  Lastly,  we  examine  the  experimental  consequences  of  our  results. 


1.  LINEAR  THEORY 

The  wave  equation  for  simple  harmonic  waves  p(x,y,z)e'“’'  propagating  in  an 
acoustic  duct  is 


^0x^  5y^ 


52  A 


az" 


p  + 


(4.1.1) 


where  p  is  the  pressure,  c  is  the  thermodynamic  speed  of  sound  in  the  medium, 
and  (0  is  the  angular  frequency  of  the  wave.  Setting  the  z  axis  parallel  to  the  axis 
of  the  duct  we  seek  a  solution  of  the  form 
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p  =  T(x,y)e''^  . 


(4.1.2) 


where  the  factor  T  satisfies  the  equation 


+ 


2  A 


ax  ay 


4'  +  7^T  =  0  , 


with 


(4.1.3) 


(4.1.4) 


the  waveguide  dispersion  relation  for  o  (k). 

For  a  rigid  wall  duct,  the  fluid  velocity  normal  to  the  surface  wall  must  be  zero. 
From  the  linearized  momentum  equation 


a  harmonic  wave  satisfies 

ipcou  =  Vp  (4.1.6) 

where  p  is  the  density  and  u  is  the  fluid  velocity.  Thus,  in  terms  of  the  acoustic 
pressure  the  appropriate  boundary  condition  for  a  rigid  wall  is 


(4.1.7) 


for  every  point  on  the  wall  perimeter  with  n  the  normal  to  the  wall. 
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There  exists  a  discrete  set  of  eigenvalues  for  which  Eqs.  (4.1.3)  and  (4.1.7) 
are  satisfied.  From  Eq.  (4.1.4)  it  can  be  seen  that  to  each  mode  corresponds  a 
specific  cutoff  frequency  cy  j  .  For  frequencies  co  greater  than  cyj,  k  is  a  real 

number  and  the  wave  propagates  along  the  tube.  For  frequencies  less  than  cyj, 
k  is  imaginary  and  the  mode  does  not  propagate. 

The  lowest  eigenvalue  y  o  corresponds  to  the  fundamental  plane  wave  mode. 
For  this  mode,  the  phase  velocity  is  the  speed  of  sound  in  the  medium.  At 
frequencies  below  the  cutoff  frequency  of  the  first  higher  mode  only  the 
fundamental  plane  wave  mode  propagates  without  attenuation. 

As  the  frequency  is  increased  above  the  cutoff  frequency  of  the  first  higher 
mode,  the  plane  wave  mode  continues  to  propagate  together  with  the  first  higher 
mode. 

The  phase  velocity  of  the  jth  mode  is 


(4.1.9) 


and  the  group  velocity  for  a  group  of  waves  centered  about  frequency  w  is 


V 


g 


5co 


(4.1.10) 


This  velocity  is  less  than  c  for  frequencies  above  cutoff. 
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For  a  waveguide  of  circular  cross  section  of  radius  R,  the  transverse 
eigenfunctions  and  eigenvalues  of  Eq.  (4.1.3)  are  (Morse  and  Ingard,  1968) 


vp  = 
m.n 


(4.1.11) 


_  7IW 

Y  m,n  ~  p  ' 


(4.1.12) 


where  r  and  (j)  are  the  polar  coordinates  about  the  central  axis  of  the  cylinder 
and  where  each  pair  of  subscripts  m,  n  defines  a  wave  mode  propagating  in  the 
waveguide.  The  subscript  m  indicates  the  number  of  nodal  planes  while  the 
subscript  n  indicates  the  number  of  nodal  cylinders.  The  alternative  possibilities 
of  using  cosine  or  sine  functions  express  the  fact  that  the  non  axially  symmetric 
modes  are  degenerate. 

For  each  integer  value  of  m  there  exist  a  set  of  values  w  which  satisfy  the 
boundary  condition 

Ji,(TO)  =  0.  (4.1.13) 

The  solutions  w  which  satisfy  this  condition  are  denoted  as  w  =  The  plane 
wave  mode  corresponds  to  ago  =0.  Jo(^«o,o)  =  I-  f^^om  Eq.  (4.1.12),  the  cutoff 
frequency  of  the  (m,n)th  mode  is  afn.nC/2R  Hz. 

The  first  mode  above  the  plane  wave  mode  is  the  (1 ,0)  mode.  This  mode  has 
one  nodal  plane  extending  radially  out  from  the  axis  and  no  nodal  circles.  Since 
the  axial  velocity  distribution  is  proportional  to  Ji(7iai  q/  R)cos((t)),  the  magnitude 
of  the  axial  velocity  along  a  line  transverse  to  the  nodal  plane  is  0  at  the  node 
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(  (|>=±7i/2)  and  a  maximum  at  r=R.  However,  the  axial  velocities  in  the  two  halves 
of  the  waveguide  are  opposite  in  phase.  This  dipole  mode  is  the  dominant  mode 
excited  when  the  cylindrical  waveguide  is  driven  above  cutoff  with  two  speakers 
in  an  anti-symmetric  push-pull  configuration. 

Generally  the  wavenumber  k  may  have  real  and  imaginary  parts  and  we  can 
write 

k  =  ^  +  i7rTij  ,  (4.1.14) 

where  Xj  is  the  wavelength  along  the  axis  of  the  waveguide  of  the  jth  mode  and 
Tij  is  an  attenuation  factor.  The  sources  of  the  attenuation  are  viscous  and 
thermal  losses  at  the  duct  walls  as  well  as  losses  due  to  the  acoustic  admittance 
of  the  boundary  surface.  For  a  plane  wave  mode  the  attenuation  due  to 
thermoviscous  losses  is  (Landau  and  Lifshits,  1959) 

no,o  =  ;^(l  +  VS^U-l))  .  (4.1.16) 

where  v  is  the  kinematic  viscosity  coefficient,  Pr  is  the  PrantI  number,  and  ^  is 
the  ratio  of  specific  heats.  An  expression  for  the  attenuation  of  the  higher  order 
propagating  mode  (1,0)  as  measured  along  the  z  axis  is 


hi.o 


^/^Rc 


:[l-(l/:rta,o)^ 


+VpF(5-i) 

> 


(4.1.16) 
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The  modification  of  the  term  above  from  its  form  in  (4.1.15)  can  be  understood 
by  recognizing  that  the  distance  traveled  by  the  wavefronts  in  the  direction  of 
their  wavevectors  is  (Cp/c)  times  greater  than  the  distance  measured  along  the  z 
axis.  Because  of  the  circular  geometry  the  wavefronts  are  not  straight.  This 
effect  is  .  .^presented  by  the  factor  .  Thus,  the  viscous  losses  are 

multiplied  by  (Cp/c[1-(1/7raio)^]).  The  contribution  due  to  thermal  losses  remains 
the  same  because  they  are  independent  of  the  shear  component  of  the  velocity. 

The  distance  at  which  the  pressure  field  decreases  by  a  factor  e”^ ,  the  e 
folding  distance,  is  the  inverse  of  the  expression  (4.1.16).  It  approaches  0  as  the 
frequency  approaches  the  mode's  cutoff  frequency,  it  has  a  maximum  when 
(0  /  cyi  0  is  approximately  equal  to  3’'^,  and  goes  as  for  high  frequencies. 

2.  THE  NONLINEAR  COEFFICIENT 

The  combined  effects  of  nonlinearity  and  dispersion  may  force  a  single 
wavetrain  to  split  into  two  separate  wavetrains  moving  with  unequal  velocities.  A 
calculation  of  the  magnitude  of  the  nonlinear  coefficient  is  thus  required  to 
determine  whether  the  effect  is  measurable  for  the  specific  case  of  the  (1,0) 
cylindrical  waveguide  mode. 

A  general  expression  for  the  nonlinear  coefficient  for  the  wavenumber  kj  has 
been  obtained  by  Keller  and  Millman  (1970).  They  consider  the  effect  of  a  small 
finite  amplitude  sound  wave  propagating  in  a  rigid  cylindrical  waveguide.  A 
perturbation  method  is  used  by  expanding  the  velocity  potential  O  and  the 
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pressure  p  in  powers  of  the  dimensionless  amplitude  e  to  solve  the  equations  of 
fluid  dynamics. 

Specifically,  consider  the  Bernoulli  equation 


1 

- +  — 

8t  2 


dp 

p(p) 


=  0. 


(4.2.1) 


and  the  continuity  equation 


Pp[p(P)K<l>-^  =  2V<t..yf^]+i(V®  V)(V(t>f  ,  (4.2.2) 

dl  \  cl  J  ^ 

where  O  is  the  velocity  potential,  p  is  the  pressure  p(p)  is  the  density,  p^°*  is  a 
constant  and  Pp  is  the  derivative  of  the  pressure  with  respect  to  the  density. 

Equations  (4.2.1)  and  (4.2.2)  must  be  supplemented  with  the  boundary 
condition 


(4.2.3) 


corresponding  to  the  vanishing  of  the  normal  component  of  the  velocity  at  the 
waveguide  perimeter. 

Regarding  the  frequency  <o  as  given  and  k  as  unknown  and  taking  the  z  axis 
as  parallel  to  the  waveguide  axis,  solutions  of  the  form 

<P  =  ci)(x,y,kz-cot),  p  =  p(x,y,kz-(Dt),  (4.2.4) 


are  sought. 
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Introducing  the  variable  change  z'  =  kz  and  t'  =  wt  and  then  dropping  the 
primes  from  the  resulting  equations  allows  us  to  rewrite  Eqs.  (4.2.1)  and  (4.2.2) 
in  the  form 


Cl) 


dt 


1 

+  — 

2 


j 


,az  j 


dp 

P(P) 


=  0, 


(4.2.5) 


Pp[p(P)] 


ax  ay' 


+  k' 


a^cp 

dz^ 


-CO 


2  a^o 


1  a  ao  a  acp  a 

2 a  ax  ax  ay  ay 


az  az 


^a^cp 


a^cp  . 

+ - r-  +  k 


ax^  ay' 


2  a^cp 

az^ 


where 

cp(x,y,z-t  +  27i)  =  0(x,y.z-t), 
p(x,y,z-t  +  27i)  =  p(x,y,z-t). 


(4.2.6) 


(4.2.7) 


The  Taylor  expansion  of  <P,  p.  and  k  in  terms  of  e  can  be  obtained  as  follows 
cp„(x,y,z-t,E)  =  ECpJ,^^+y(pJ,2^+...  ,  (4.2.8) 

p„(x,y,z-t,E)  =  p^°^+epI,^^  +  yplf^+...  ,  (4.2.9) 


k„(e)  =  k'”'+ek<«+^k®+...  , 


(4.2.10) 


where  the  superscripts  denote  the  jth  derivative  with  respect  to  e  and  the 
subscript  denotes  the  propagating  mode.  The  justification  for  Eq.  (4.2.10)  is  that 
the  Taylor  expansions  of  O  and  p  in  terms  of  e  do  not  exist  unless  k  also 
depends  on  e  . 
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Keller  and  Millman  (1970)  insert  Eqs.  (4.2.8)  -  (4.2.10)  into  Eqs.  (4.2.5)  and 
(4.2.6).  They  show  that  =  0  and  obtain  the  following  general  expression  for 
the  coefficient 

2o)|[l<f  [  4.2h|p„4'„  +q„'C=)')'„d*dy  .  (4.2.11) 

where  c  is  the  speed  of  sound  in  the  medium,  k|,°*  is  the  linear  wavenumber  for 

2  2 

which  k[,°-  =  ^-Yn  (  Yn  is  the  eigenvalue  of  the  wave  equation  (4.1.3)  and 

c 

Tp  is  the  corresponding  eigenfunction),  is  the  transverse  Laplacian,  is 
determined  by  solutions  to 

+4y;)p„  +2hj'1'„4  ,  (4.2.12) 

and 


For  an  ideal  gas  equation  of  state  p(p)=Ap' 
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(4.2.14) 


r 


and 


(4.2.15) 


We  now  specialize  these  results  to  the  case  of  the  first  non-axially  symmetric 
mode  propagating  in  a  cylindrical  waveguide  of  radius  R.  From  Eq.  (4.1.11)  the 
normalized  solution  of  the  full  wave  equation  for  this  mode  is 


V^a^o  cos(<t))J^ 


^1,0  = 


Ttaior 

R 


RJi(7taio)^(7raio)' 


(4.2.16) 


This  mode  is  doubly-degenerate  and  even  when  an  external  drive  breaks  the 
symmetry,  nonlinearities  resonantly  couple  the  two  eigenmodes.  However, 
because  a  real  waveguide  is  not  perfectly  circular  the  degeneracy  is  removed 
and  the  frequency  difference  is  a  monotonic  function  of  the  eccentricity.  Thus 
we  will  assume  a  nondegenerate  solution  of  the  form  (4.2.16)  and  determine  the 
nonlinear  coefficient  (4.2.1 1)  accordingly. 

To  obtain  an  expression  for  q  insert  (4.2.16)  into  (4.2.12)  and  transform 
to  polar  coordinates.  The  right  hand  side  of  the  inhomogeneous  Eq.  (4.2.12)  is 
of  the  form  F(r)  +  H(r)  cos(2<j)),  where  F(r)  and  H(r)  are  in  terms  of  products  of 
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two  Bessel  functions.  The  solution  of  (4.2.12)  satisfying  the  boundary  condition 
(4.2.3),  which  becomes  P,  o^(R)  =  0,  is 


R  R 

p,o(r)  =  jGo(r,s)F(s)sds  +cos(2<t))jG2(r,s)H(s)sds  .  (4.2.17) 

0  0 

Here,  the  Green's  functions  Gm(r,s)  are  given  by 

TT  Jm(2na|.or< /R) 

2  J„’(27ta,o)  (4.2.18) 

X  [j^(27ta,or> /R)N^’(2m,o)-N^(27ta,or^ /R)Jn,'(27ia,o)]  . 

where  is  the  Neumann  function  of  order  m,  and  r>  (r<)  stands  for  the  greater 
(smaller)  of  r  and  s.  Using  (4.2.14)-(4.2.18)  we  obtain  for  kjJJ  the  following 


kS(<o)  = 


■J^a 


1,0 


RJ|(^ai,o)-/(  7ta,y)^  -1 

{o.305  +  0.38x-  -1.463X''  +0.371x^} 


(4.2.19) 


7t^a,o‘ 


2c2R'kS'3 


whei'e  X  =  C(>R/(C7rai  g)  is  the  ratio  of  the  frequency  to  the  cutoff  frequency.  The 
expression  (4.2.19)  has  meaning  for  x>1.  It  is  negative  for  1<  x  <  1.897  and 
becomes  positive  for  frequencies  larger  than  1.897  times  bigger  than  the  cutoff 
frequency. 
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The  wavenumber  is  the  result  of  developments  in  powers  of  e  which  is  the 
amplitude  of  the  linearized  solution.  In  terms  of  the  ratio  of  the  acoustic  pressure 
to  the  ambient  pressure  we  obtain 


'(0‘  (7ra,o)^' 

c‘ 

'RJ,(m|o)-y/(7ta,o)^  -1  ’ 

u(2) 

'^1,0 

fp] 

2 

kPo  > 

The  propagation  constant  q  depends  upon  amplitude.  It  is  less  than  the  linear 
value  for  frequencies  between  the  cutoff  frequency  and  1 .897  times  the  cutoff 
frequency.  According  to  the  physical  arguments  presented  in  the  Chapter  1 ,  a 
modulation  of  a  wavetrain  will  be  stable  and  an  intial  disturbance  will  split  into 
two  disturbances  moving  with  different  group  velocities.  On  the  other  hand,  for 
frequencies  greater  that  1 .897  times  the  cutoff  frequency  the  phase  velocity  is  a 
decreasing  function  of  amplitude  and  modulations  of  a  wavetrain  will  become 
unstable  to  the  formation  of  solitons. 

The  formalism  developed  in  the  previous  chapters  relies  on  the  determination 
of  the  nonlinear  coefficient  ©2.  Equation  (4.2.20)  gives  the  fequency 
dependence  of  the  propagation  constant  k,  q  .  Thus,  we  can  solve  implicitly  for 
the  frequency  as  a  function  of  wavenumber  k  =  k^  q  •  If  we  set  co=a)o+o)2(p/po)^,  it 
can  be  shown  that 


(o,(k)  =  -<oUk)- 


RJi(mio)^/(  7ra,Q)^  - 1 
V?7ra,o^c0o(k) 


kS(cOo(k)) 


Ipo> 


(4.2.21) 
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where 


coo(k)  =  c 


Tia 


1,0 


+k' 


(4.2.22) 


is  the  linear  dispersion  relation  and  (Hq  is  the  linear  group  velocity  at  the 
wavenumber  k. 


3.  EXPERIMENTAL  CONSEQUENCES 

For  a  fixed  frequency,  Eq.  (4.2.20)  predicts  that  the  wavelength  of  a  nonlinear 
wave  changes  with  amplitude.  It  increases  with  amplitude  for  frequencies  less 
than  1.897  times  the  cutoff  frequency  and  decreases  with  amplitude  for 
frequencies  that  are  larger  than  1.897  times  the  cutoff  frequency.  The  mode  is 
modulationally  stable  in  the  former  case  and  we  will  restrict  our  analysis  to  this 
frequency  range. 

The  effect  of  increasing  the  wavelength  with  amplitude  can  be  determined 
experimentally  by  measuring  the  relative  phase  between  two  points.  If  the 
distance  d  between  two  points  is  less  that  the  wavelength 


Acp  _  Ak 
cp  k 


(4.3.1) 
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where  A(p  is  the  change  in  phase  difference  (p  between  two  points  due  to  a 
change  in  amplitude,  and  Ak  =  k,o  -kjo-  The  effect  can  be  magnified  if  the 

distance  between  the  two  points  several  wavelenghts  apart.  In  this  case 


Acp  _  Ak 
9  +  27iN  k 


(4.3.2) 


where  N  is  the  integer  number  of  wavelengths  in  the  distance  d  between  the  two 
points.  The  slope  of  the  measured  relative  phase  between  two  points  as 
function  c‘  the  square  of  the  amplitude  will  thus  yield  an  experimental  value  of 
the  nonlinear  coefficient.  This  value  in  turn,  can  be  used  to  determine  the 
separation  in  time  as  a  function  of  distance  of  the  split  pulses. 

To  determine  the  difference  of  the  times  of  arrival  at  a  given  location  of  the 
two  disturbances  that  result  from  the  splitting  of  an  initial  disturbance  we  can  use 
Eq.  (3. 1 .4).  In  terms  of  the  ratio  p/po  of  the  acoustic  pressure  to  the  ambient 
pressure,  (3.1.4)  becomes 


^  =  o,„'(k)±{a),(k)c»o"(k)}''’ 


(4,3.3) 


A  finite  extent  modulation  of  a  nonlinear  wavetrain  generated  at  a  source  at  z=0 
will  split  into  two  disturbances  that  arrive  at  a  distance  z  with  a  time  difference 
given  by 


Atn(Z)  =  Z 


2(coo  CO2) 


U2 


P_ 

VPo 


(4.3.4) 
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In  order  to  determine  if  this  effect  is  observable,  we  must  also  determine  if  the 
time  spread  due  to  dispersion  is  smaller  than  the  time  due  to  splitting. 

The  average  velocity  of  a  smooth  pulse  or  a  finite  wavetrain  that  does  not  cut 
off  too  violently  is  the  group  velocity  Vg=dcoo/dk=wo^  Due  to  dispersion  the  pulse 
spreads  because  different  frequency  components  move  with  different  phase 
velocities.  The  spreading  of  the  pulse  can  be  accounted  for  by  noting  that  a 
pulse  with  initial  temporal  width  Ato  must  have  an  inherent  spread  of  frequencies 
Aco  ~  (1/Ato).  This  means  that  the  group  velocity  spread  is  of  the  order 


Av, 


Wo 


Aw  - 


WgAto 


(4.3.5) 


which  at  a  point  z  implies  a  spiead  in  time  of  the  order  of  AVg  z  /  (cocY.  Because 
the  uncertainty  in  time  is  the  square  root  of  the  sum  of  squares  the  width  At(z)  at 
location  z  is 


Atd(z) 


(Ato)' 


Wn  Z 


W 


0  Ato  / 


(4.3.6) 


The  expression  for  At(z)  shows  that  in  a  medium  with  dispersion  a  narrow  pulse 

spreads  rapidly  because  of  its  broad  spectrum  in  frequencies.  Thus  the  criterion 

2 

for  a  small  change  in  shape  is  that  Wq  Ato  »  • 

Thus,  for  the  splitting  to  be  observable  fora//  distances  we  require  that 
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2  AtgCOo 


f  CO,  1" 


[<^o"  ) 


(4.3.7) 


For  a  10  ms  disturbance  superimposed  on  a  signal  with  a  frequency  10%  higher 
than  the  cutoff  frequency  the  minimum  ratio  of  acoustic  pressure  to  ambient 
pressure  is  0.005,  or  148  dB  (re  20  pPa). 

The  estimate  (4.3.7)  guaranties  that  the  curves  (4.3.4)  and  (4.3.6)  do  not 
intersect.  For  low  intensity  levels  the  two  curves  intersect.  However,  it  should 
still  be  possible  to  observe  the  splitting  of  a  pulse  into  two  pulses  for  distances 
less  than  the  distance  given  by  the  real  solution  to 


(Ato)'  + 
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(4.3.8) 


provided  the  value  At^  given  by  (4.3.4)  is  within  the  limits  of  instrumental 
resolution.  Spreading  due  to  dispersion  will  dominate  for  distances  larger  than 
the  distance  determined  by  the  real  solution  of  (4.3.6).  However,  even  when  the 
splitting  of  a  disturbance  into  two  cannot  be  observed,  the  effects  of  the  double 
group  velocity  will  be  to  spread  a  given  pulse  beyond  the  values  predicted  by 
linear  dispersion. 
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CHAPTER  5.  EXPERIMENTAL  RESULTS 


In  this  chapter  we  discuss  the  experimental  results  related  to  propagation  of 
the  (1 ,0)  cylindrical  waveguide  mode  in  the  frequency  regime  for  which  the  mode 
is  modulationally  stable.  We  briefly  describe  the  experimental  apparatus  used  in 
our  experimental  investigations.  The  experimental  results  include 
measurements  of  the  (1 ,0)  mode  linear  cut  off  frequency,  group  velocity,  and 
attenuation  constants  for  frequencies  above  the  cutoff  frequency.  From 
frequency  response  measurements  we  also  determine  the  relative  amplitudes  of 
the  plane  wave  and  dipole  modes.  Also,  we  measure  the  spreading  of  a  low 
amplitude  acoustic  pulse.  Finally,  investigations  of  nonlinear  effects  including 
the  determination  of  the  nonlinear  coefficient  at  4200  Hz  and  preliminary  time 
domain  observations  of  group  velocity  splitting  in  amplitude  modulated  and  in 
pulse  modulated  acoustic  waves  are  presented. 

1.  EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  is  a  traveling  wave  tube  made  of  an  extruded 
aluminum  pipe  of  average  inner  diameter  1.957  inches,  outer  diameter  2.377 
inches,  length  70  feet,  which  consist  of  seven  10  foot  sections  joined  by  collars 
and  flanges.  An  anechoic  end  termination  made  of  steel  wool  with  a  tapered 
density  distribution  spans  the  length  of  the  7th  section.  Absorber  and  junction 
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reflection  amplitudes  are  typically  -40dB  down  for  frequencies  between  1  kHz  and 
4kHz  (Dorff,  1991). 

Each  of  the  first  six  sections  of  pipe  have  4  microphone  ports  which 
accommodate  ENDEVCO  series  851  OB  piezoresistive  pressure  transducers. 

The  ports  have  a  uniform  spacing  of  2.5  feet.  An  additional  port  was  added 
directly  below  the  2nd  port  in  the  5th  section  (from  the  drivers'  end  of  the  pipe). 
This  port  adds  the  capability  of  instrumentally  eliminating  the  detection  of  the 
plane  wave  mode  by  connecting  the  microphones  with  opposite  polarity  and 
adding  their  output. 

Two  JBL  model  2445H  compression  drivers  were  mounted  at  one  end  of  the 
wave  tube  via  a  two-to-one  adapter.  The  first  microphone  port  is  at  a  distance  of 
0.83  m  from  the  drivers.  The  drivers  are  titanium  diaphragm  drivers  rated  at  100 
watts  from  500  Hz  to  1000  Hz  and  150  watts  above  1  kHz. 

A  two-to-one  adapter  was  designed  such  that  the  drivers  could  be  mounted 
directly  opposed  (Fig.  5.1).  The  inner  diameter  of  the  adapter  was  machined  to  a 
diameter  of  1 .975  inches  to  match  the  inner  end  diameter  of  the  wave  tube.  The 
speakers  were  mounted  flat  to  the  adapter  onto  rectangular  attachment  flanges. 
A  2  inch  to  1  inch  tapered  bore  made  through  the  attachment  flanges  and  the 
adapter  connected  the  2  inch  driver  opening  to  the  1.975  inch  diameter 
waveguide  opening. 

In  this  configuration,  the  drivers  operate  in  push  -  pull  when  driven  180 
degrees  out  of  phase.  The  plane  wave  mode  can  in  principle  be  nulled  by 
adjusting  the  amplified  signal  to  each  driver.  In  this  manner,  at  frequencies 
above  cutoff  the  (1,0)  mode  is  the  main  propagating  mode. 
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2.  EXPERIMENTAL  RESULTS 

a)  Frequency  Response 

In  an  acoustic  waveguide  the  plane  wave  mode  propagates  at  all 
frequencies.  Although  the  experimental  procedure  calls  for  nulling  the  plane 
wave  component  of  the  acoustic  field,  the  frequency  responses  of  the  drivers  are 
not  identical.  Therefore,  even  if  the  plane  wave  component  is  nulled  at  one 
frequency,  it  may  not  be  so  for  other  frequencies  in  the  sweep. 

The  role  played  by  the  plane  wave  field  was  investigated  by  taking 
frequency  response  measurements  with  a  swept  sine  source  at  the  first  port  0.83 
m  from  the  drivers'  end  and  at  each  alternate  port.  The  frequency  span  of  the 
sweep  is  from  3900  to  4400  Hz,  with  a  resolution  of  401  points  per  sweep.  The 
plane  wave  component  was  nulled  at  3900  Hz. 

To  investigate  the  effects  that  the  presence  of  a  plane  wave  might  have  in 
the  dipole  field,  we  obtained  frequency  responses  measured  at  the  double  ports 
located  at  13.79  m  from  the  drivers'  end.  The  double  port  arrangement  allows 
measurements  to  be  obtained  by  the  output  of  a  single  microphone  as  well  as 
the  output  of  the  sum  of  the  two  microphones  with  the  same  and  with  opposite 
polarities. 

The  frequency  response  measured  from  a  single  microphone  is  the 
superposition  of  a  plane  wave  and  a  dipole  field.  For  a  given  frequency,  the 
corresponding  wavelengths  of  the  two  modes  are  not  the  same.  Consequently, 
the  relative  phase  of  the  corresponding  acoustic  pressure  is  a  function  of 
frequency.  Thus  we  should  expect  that  the  interference  between  the  two  fields 
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at  a  given  point  be  reflected  in  the  frequency  response  at  this  point.  Fig.  5.2a 
shows  a  sinusoidal  variation  of  the  amplitude  as  a  function  of  frequency. 

Fig.  5.2b  shows  the  output  at  the  dual  port  of  two  microphones  with  the 
same  polarity.  In  this  arrangement  the  dipole  field  is  canceled  and  the  output 
corresponds  to  a  frequency  response  of  the  plane  wave.  The  inconsistency  in 
the  amplitudes  for  the  frequency  responses  of  Fig.  5. a  as  compared  to  that  in 
Figs.  5.b,c  is  due  to  the  difference  in  the  microphone  gain  used  to  obtain  these 
signals.  The  microphone  gain  in  Fig.  5. a  is  500  while  the  gain  used  to  obtain 
Figs.  5.b,c  is  10k.  Note  that  for  frequencies  above  4050  Hz,  the  amplitude 
decreases  with  increasing  frequency.  In  Figs.  5.2a,b  the  frequency  response 
abruptly  dips  at  4030.2  Hz  which  corresponds  to  the  cutoff  frequency  for  the 
(1,0)  mode. 

Alternatively,  Fig.  5.2c  shows  the  output  at  the  dual  port  of  two 
microphones  with  opposite  polarity.  In  this  arrangement  the  plane  wave  field  is 
canceled  and  the  output  corresponds  to  the  frequency  response  of  the  (1 ,0) 
mode  dipole  field.  For  frequencies  greater  than  4200  Hz,  the  frequency 
response  is  essentially  flat. 

A  comparison  of  Figs.  5.2a,b&c  leads  us  to  the  conclusion  that 
interference  is  occurring  between  the  plane  wave  and  the  (1,0)  mode.  As  shown 
in  Figs.  5.2b&c  the  plane  wave  mode  is  observed  to  be  between  9  and  25  dB 
down  from  the  dipole  mode  for  frequencies  above  the  cutoff. 


Amplitude  (dBVrmt)  Amplitude  (dBVrms) 


•30 


Frequency  (Hi) 


Frtqucncy  (Hi) 


Frequtney  (Hi) 


Fig.  5.2  Frequency  response  of  the  traveling  wave  tube  measured  at 
the  dual  port  13.79  m  from  the  drivers,  a)  Output  from 
single  microphone.  Mike  gain  is  500.  b)  Output  from  two 
microphones  with  equal  polarity  c)  Output  from  two 
microphones  with  opposite  polarity.  Mike  gain  for  b,c  is  10k. 
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b)  Attenuation  Constant 

From  the  frequency  response  at  each  port  we  obtained  microphone  signal 
amplitudes  for  401  frequencies  between  3900  Hz  and  4400  Hz.  A  matrix  was 
constructed  from  this  data.  The  rows  of  the  matrix  contain  the  signal  amplitudes 
at  fixed  frequency  for  14  microphone  port  locations.  The  columns  list  signal 
amplitudes  at  a  fixed  location  for  401  separate  frequencies  between  3900  and 
4400  Hz.  Using  this  matrix,  we  obtained  plots  of  the  amplitude  versus  distance 
for  each  frequency.  From  an  exponential  fit  performed  on  the  data,  we  obtained 
the  attenuation  constants  for  frequencies  between  4050  and  4400  Hz. 

Fig.  5.3a  is  a  plot  of  the  e'^  folding  distance  as  a  function  of  frequency. 
The  attenuation  constant  corresponds  to  the  inverse  of  this  distance.  The  plot  of 
this  data  shows  irregularly  high  and  low  e'^  folding  distances  between  the 
frequencies  4100  to  4200  Hz.  As  suggested  by  Fig.  5.2  the  variation  in  the  e'^ 
folding  distance  between  these  frequencies  is  due  to  the  rapid  variations  in  the 
frequency  response  at  the  different  ports  between  these  frequencies.  That  is, 
the  frequency  corresponding  to  a  peak  at  one  port  may  not  necessarily 
correspond  to  a  peak  at  another  port.  This  observation  is  supported  by  the  fact 
that  variation  in  the  e'^  folding  distances  is  not  observed  at  frequencies  where 
the  amplitude  of  the  frequency  response  is  relatively  flat. 

In  Fig.  5.3b  we  have  removed  the  extreme  high  and  low  attenuation 
factors  in  the  frequency  range  4100  to  4200  Hz.  The  solid  curve  in  Fig.  5.3b, 
obtained  from  Eq.  (4.1.15),  corresponds  to  the  predicted  attenuation  length  due 
to  thermoviscous  losses  .  Excellent  agreement  between  theory  and  experiment 
is  observed  up  to  4300  Hz  where  pipe  radiation  losses  possibly  become 
important. 
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Fig.  5.3  Plot  of  the  1/e  folding  distance  as  a  function  of  frequency. 
The  extreme  high  and  low  variations  in  (a)  have  been 
removed  in  (b).  The  solid  line  corresponds  to  theory. 
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c)  Group  Velocity  Measurements 

The  frequency  response  plots  of  Fig.  5.2  show  a  change  in  behavior  for 
frequencies  above  the  cutoff  frequency  of  4030.2  Hz.  According  to  the  theory, 
when  the  drivers  are  in  a  push-pull  configuration  the  dominant  propagating  mode 
is  the  (1 ,0)  mode.  For  this  mode,  the  group  velocity  depends  on  frequency 
according  to  the  relation  (4.1.10) 

Vg=c(l-(co,/co)')'''  (5.2.1) 

where  is  the  cut  off  frequency  of  the  (1 ,0)  mode.  Thus  for  a  fixed  distance  D, 
the  time  of  travel  At  is 

(l/At)^=cVD={l-a),=  (l/«)")}  .  (5.2.2) 

We  made  measurements  of  the  time  of  flight  of  a  100  cycle  burst  with 
carrier  frequencies  above  cutoff.  The  envelope  was  sinusoidal.  The  time  delay 
measurements  were  made  at  a  single  port  and  a  fixed  location  13.79  m  from  the 
drivers  on  a  Nicolet  Oscilloscope  triggered  externally  from  the  driver  signal. 

Time  delays  were  measured  from  the  time  of  trigger  to  the  time  of  half  maximum 
signal  amplitude  for  five  frequencies  (Coppens,  1966).  Fig.  5.4  shows  the  plot  of 
the  resulting  data.  As  predicted  from  (5.2.2)  there  is  a  linear  relationship 
between  the  inverse  of  the  square  of  the  time  delay  ( 1/At^ )  and  the  inverse  of 
the  square  of  the  frequency  ( 1/f^ ).  A  least  square  fit  of  the  data  yields  a  cutoff 
frequency  of  4040  Hz.  This  result  is  within  0.2%  of  the  value  obtained  from  the 
frequency  response  data.  The  experimentally  obtained  value  for  the  speed  of 
sound  is  338.7  m/s. 
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Fig.  5.4  Group  velocity  measurements.  The  graph  shows  a  linear 
dependence  between  the  inverse  of  the  square  of  the  time 
delay  and  the  inverse  of  the  square  of  the  frequency  for  the 
(1 ,0)  mode  for  5  frequencies  above  cutoff.  Slope  =  877.9  ; 
Intercept  =  -1.43  x  10^°. 
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d)  Measurements  of  Pulse  Spreading 

Equation  (4.3.6)  estimates  the  spreading  of  an  acoustic  pulse  due  to 
dispersion.  This  equation  predicts  that  the  spreading  is  significant  only  at 
distances  much  greater  than  the  length  of  our  waveguide  or  for  frequencies  very 
close  to  cutoff.  Specifically,  predictions  for  the  spreading  of  a  21 .7  ms  pulse 
propagating  a  distance  of  7.69  m  range  from  0.002  ms  for  a  4.6  kHz  signal  to  0.1 
ms  for  a  4.2  kHz  signal.  These  times  are  within  the  instrumental  uncertainty  of 
the  Nicolet  oscilloscope. 

In  order  to  verify  the  prediction  that  the  spreading  due  to  dispersion  was 
negligible,  we  analized  the  speeding  of  a  low  amplitude  N  cycle  pulse  of  21 .7  ms 
duration.  The  carrier  was  modulated  to  obtain  a  smooth  sinusoid  envelope.  The 
number  of  cycles  N  was  varied  to  keep  the  input  pulse  duration  constant. 
Measurements  of  pulse  spreading  were  made  by  measuring  the  total  time 
duration  of  the  pulse  signals  obtained  at  locations  0.83m  and  7.69m  from  the 
drivers.  Table  5.1  shows  the  measurements  obtained  for  5  frequencies. 

Table  5.1 

4.6  kHz  4.5  kHz  4.4  kHz  4.3  kHz  4.2  kHz 


(ms) 

(ms) 

(ms) 

(ms) 

(ms) 

0.0  m 

21.5  ±  .2 

21.7  ±  .2 

21.7  ±  .2 

21.7  ±  .2 

21.7  ±  .2 

0.83  m 

22.0  ±  .4 

24.0  ±  .4 

24.7  ±  .4 

25.5  ±  .2 

27.1  ±  .2 

7.69  m 

21.6  ±  .4 

21.8  ±  .4 

22.0  ±  .4 

22,2  ±  .4 
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The  data  verify  that  for  the  propagation  distances  and  frequencies 
measured  the  contribution  to  the  spreading  from  dispersion  is  negligible. 

The  spreading  of  the  pulse  at  0.83  meters  is  apparently  due  to  the  fact 
that  at  this  distance  the  plane  wave  has  not  yet  unraveled  from  the  center  of  the 
dipole  mode  pulse.  That  is,  the  microphone  signal  consists  of  a  plane  wave 
component  as  well  as  a  dipole  component  where  each  travel  with  different 
velocities  of  propagation  and  therefore  arrive  at  the  microphone  at  different 
times.  For  example,  in  considering  the  data  for  the  4.4  kHz  signal,  the  plane 
wave  component  arrives  at  the  microphone  in  2.45  ms.  The  dipole  component 
arrives  in  6.1  ms.  The  combined  signal  would  therefore  appear  to  be  3.6  ms 
longer  than  the  input  signal.  Experimentally,  the  signal  at  the  .83  meter  location 
for  4.4  kHz  is  observed  to  be  3  ms  longer  than  the  input. 

Signal  distortion  prevented  a  reliable  measurement  of  the  duration  of  the 
4.2  kHz  signal  from  the  7.69m  location. 
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e)  Measurement  of  the  Nonlinear  Coefficient 

As  shown  in  Fig.  5.4,  for  frequencies  above  cutoff  the  mode  is  dispersive. 
For  small  but  not  infinitesimal  amplitude,  the  dispersion  relation  for  the  (1,0) 
mode  is  of  the  form  given  by  Eq.  (2.2.1).  According  to  this  expression,  for  fixed 
frequency  a  change  in  the  amplitude  induces  a  change  in  wavelength  of  the 
mode.  If  the  nonlinear  coefficient  is  positive  the  changes  in  the  wavenumber 
decrease  linearly  with  the  square  of  the  amplitude.  Thus,  by  measuring  the 
relative  phase  between  two  points  as  a  function  of  amplitude  we  can  determine 
both  the  validity  of  Eq.  (2.2.1)  and  the  nonlinear  coefficient. 

Using  a  PAR  5210  lock-in  amplifier  we  determined  the  amplitude  of  the 
acoustic  field  at  13.79  m  from  the  drivers  and  the  phase  of  this  signal  relative  to 
a  reference  signal  located  at  a  distance  of  9.98  m  from  the  drivers.  By  varying 
the  driver  amplitude  at  a  fixed  frequency,  in  this  case  4200  Hz,  we  obtained  the 
relative  phase  between  the  two  signals  as  a  function  of  the  amplitude  squared. 
The  plane  wave  component  of  the  acoustic  field  was  nulled  at  this  frequency. 

The  data  in  Fig.  5.5  are  in  terms  of  the  acoustic  pressure  at  the  source.  A 
reasonable  estimate  of  the  amplitude  at  the  source  was  obtained  by  multiplying 
the  calibrated  acoustic  amplitude  at  13.79  m  by  the  inverse  attenuation  factor 
gi3  79/3  39  yvhere  3.39  m  is  the  1/e  folding  distance  of  a  (1,0)  mode  of  frequency 
equal  to  4200  Hz. 

Fig.  5.5  shows  a  linear  relation  between  the  relative  phase  and  the  square 
of  the  acoustic  amplitude  at  the  source.  In  order  to  get  an  idea  of  the  intensity 
levels  used  for  this  measurement,  the  pressure  amplitude  has  been  normalized 
to  the  atmospheric  pressure.  The  plotted  range  varies  between  150  dB  and  170 
dB  (re  20  pPa).  The  nonlinear  coefficient  can  be  determined  from  using  the 
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slope  of  the  line  and  Eq.  (4.3.2).  For  a  frequency  of  4200  Hz,  the  theoretical 
value  predicted  by  Eq.  (4.2.20)  yields  Ak  =  -12.58  (p/Po)  ■  In  the  span  of  3.81  m, 
corresponding  to  the  distance  between  the  two  ports  there  are  13  wavelegths. 
Thus  we  obtain 

—  ((p  +  27rN)  =  -47.03f-^l  (5-2.3) 

k  iPo 

The  predicted  slope  is  12%  bigger  than  the  experimental  value  of  -  41 .36. 
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phase  (rad) 


Fig.  5.5  Plot  of  the  relative  phase  versus  the  square  of  the 

normalized  acoustic  pressure  at  the  source.  Pressure 
amplitudes  are  normalized  to  1  atmosphere. 
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f)  Measurements  of  Group  Velocity  Splitting 

An  initial  disturbance  superimposed  on  a  signal  with  a  frequency  10% 
higher  than  the  cutoff  frequency  splits  into  two  and  the  differential  time  of  arrival 
is  approximately  30  z(p/po)  ms  at  a  distance  z  from  the  source  (Eq.  (4.3.4)). 

Fig.  5.6  shows  the  microphone  output  on  a  Nicolet  digital  oscilloscope  for  three 
different  locations  and  three  different  intensity  levels. 

The  input  signal  is  generated  using  2  HP  3314  function  generators.  A 
4500  Hz  carrier  is  amplitude  modulated  by  a  100  Hz  inverted  half-sine  wave  at  1 
sec  intervals.  The  first,  second  and  third  rows  in  Fig.  5.6  correspond  to  intensity 
levels  of  140  dB.  144  dB,  and  156  dB  at  the  source,  respectively.  At  140  dB  no 
apparent  splitting  of  the  negative  pulse  is  observed  up  to  the  shown  distance  of 
7.69,  while  the  spreading  of  the  pulse  is  noticeable.  However  at  144  dB  and  156 
dB  the  pulse  undergoes  spreading  and  a  splitting  becomes  readily  apparent  for 
distances  of  4.64  m  and  beyond. 

We  performed  quantitative  analysis  of  the  signals  displayed  in  Fig.  5.6 
The  procedure  used  consists  of  multiplying  the  data  by  the  cosine  of  the  carrier 
frequency,  taking  the  Fast  Fourier  Transform  (FFT)  of  the  product  and  then 
filtering  the  modulus  of  the  FFT  to  remove  the  high  frequency  components  of  the 
frequency  spectrum.  The  inverse  FFT  of  the  filtered  frequency  spectrum  is  then 
taken  to  recover  the  signal  modulation  exclusively. 

With  this  procedure  we  can  determine  the  time  spread  of  the  pulse  and 
the  differential  time  of  the  split  portions.  Measurements  were  made  from  the 
waveform  peaks.  For  a  sound  level  of  140  dB  a  10  ms  pulse  spreads  into  an 
1 1 .8  ms  pulse  at  4.64  m  and  into  a  14.2  ms  pulse  at  7.69  m.  Thus  the  delayed 
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.69  m 


A  10  ms  negative  modulation  at  1  sec  intervals  of  a  4.5  kHz  carrier  measured  at  different 
locations.  The  first,  second,  and  third  rows  correspond  to  140  dB,  144  dB,  and  156  dB  at 
the  source  respectiively. 


time  is  2.4  ms  and  it  is  not  at  all  described  by  dispersion.  Furthermore,  at  10.73 
m  the  pulse  undergoes  splitting  where  the  difference  between  the  two  peaks  is 
5.5  ms. 

Similarly,  for  an  intensity  of  144  dB  the  apparent  spread  of  a  10  ms  pulse 
at  4.64  m  is  14.7  ms,  while  at  7.69  m  the  spread  is  17.0  ms.  This  corresponds  to 
a  delay  of  2.3  ms.  The  split  pulses  at  these  locations  are  4.1  ms  and  6.1  ms 
apart  respectively.  Remarkably  there  is  a  delay  of  2  ms  in  the  split  pulses 
between  these  two  locations.  That  is,  the  splitting  the  pulse  occurs  with  a  delay 
nearly  equal  to  the  differential  spread. 

Furthermore,  a  comparison  of  the  144  dB  and  the  156  dB  signals  at  7.69 
m  shows  a  significant  increase  in  the  time  interval  between  two  peaks  of  the 
modulation.  Fig.  5.7  shows  that  for  the  lower  amplitude  signal,  the  time  interval 
between  the  two  most  negative  peaks  of  the  modulation  is  6.1  ms.  For  high 
amplitude  the  time  interval  between  these  same  two  peaks  is  8.4  ms.  The  time 
delay  is  2.3  ms  verifying  that  the  splitting  increases  as  a  function  of  amplitude. 

On  the  other  hand,  the  17  ms  spread  is  the  same  for  both  amplitudes. 

Using  a  single  HP  3314  function  generator  in  the  N  cycle  mode,  we 
generated  a  pulsed  sine  wave  of  50  cycles  with  frequency  4.2  kHz.  We  varied 
the  amplitude  of  the  driver  input  signal  and  obtained  time  domain  measurements 
of  a  microphone  output  signal  at  10.74  m  from  the  drivers.  The  envelope  of  the 
pulse  was  recovered  using  a  procedure  identical  to  that  described  above. 

Fig.  5.8  shows  the  modulations  of  a  150  dB  and  a  178  dB  bursts.  The 
time  intervals  between  the  two  marked  peaks  are  33.75  ms  and  35  56  m 
respectively.  The  time  delay  is  1 .8  ms. 
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Fig.  5.7  Modulation  envelope  of  a  4500  Hz  carrier  wave  at  7.69  m 
the  source  showing  peak  to  peak  time  intervals.  The 
intensity  level  at  the  source  is;  a)  144  dB,  and  b)  156  dB. 
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Fig.  5.8  Envelope  a  50  cycle  burst  at  a  frequency  4.2  kHz 

at  10.74  m  from  the  drivers.  The  acoustic  intensity  at 
the  source  is  a)  150  dB  and  b)  178  dB. 


CHAPTER  6.  CONCLUSIONS,  POTENTIAL  APPLICATIONS,  AND 

FUTURE  WORK 


1.  CONCLUSIONS 

This  work  reports  the  results  of  theoretical  and  experimental  investigations 
involving  the  phenomena  of  group  velocity  splitting.  We  have  presented  a 
quantitative  derivation  of  group  velocity  splitting  and  given  a  comprehensive 
description  of  the  process  whereby  a  modulation  of  finite  extent  splits  into  two 
separate  disturbances  traveling  with  two  different  velocities  of  propagation.  Also, 
we  have  predicted  that  for  modulationally  stable  nonlinear  dispersive  systems,  a 
consequence  of  group  velocity  splitting  is  to  introduce  a  frequency  modulation 
component  into  an  amplitude  modulated  signal. 

We  have  made  the  first  experimental  observations  of  the  change  in 
wavelength  with  amplitude  for  a  mode  above  cutoff  in  an  cylindrical  acoustic 
duct.  The  experimental  value  obtained  is  within  12%  of  the  theoretical 
prediction. 

To  our  knowledge,  we  have  made  the  first  experimental  observation  of  pulse 
splitting  in  a  nonlinear  dispersive  medium.  Experiments  show  that  a  single  10 
ms  duration  high  amplitude  (156  dB  )  pulse  propagating  7.69  m  down  the 
waveguide  splits  into  two  pulses  traveling  at  different  speeds.  Specifically,  a 
single  peaked  modulated  acoustic  signal  is  observed  to  split  into  two  peaks 
separated  in  time  by  4.1  ms  in  fair  agreement  with  theory.  Additionally,  since  the 
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spreading  of  the  pulse  at  these  distances  is  not  described  by  dispersion,  we 
conclude  that  the  observed  spreading  of  the  modulation  at  small  propagation 
distances  is  also  a  result  of  group  velocity  splitting. 


2.  POTENTIAL  APPLICATIONS 

The  importance  of  establishing  experimentally  the  effect  of  group  velocity 
splitting  lies  in  the  possible  application  to  fiber  optic  communication,  command, 
and  control  links.  For  example,  potential  applications  could  include  use  of  pulse 
splitting  in  fiber  optic  networks  requiring  increased  data  redundancy. 

The  linear  group  velocity  of  guided  waves  in  dielectric  fibers  is  determined  by 
a  combination  of  the  material  dispersion  and  geometrical  waveguide  dispersion. 
In  the  frequency  range  of  visible  light  and  for  a  small  difference  between  the  clad 
and  the  core  material,  the  material  dispersion  dominates  and  the  dispersion  is 
negative.  However  if  one  uses  cladding  having  much  smaller  index  of  refraction 
than  the  core  or  longer  wavelengths  a  positive  dispersion  results  because 
geometrical  dispersion  dominates  in  the  first  case  and  material  dispersion 
becomes  anomalous  in  the  second.  Under  these  conditions  pulse  splitting 
becomes  possible. 

For  optical  fibers  with  a  positive  dispersion,  another  potential  application  of 
the  effects  presented  in  this  work  is  the  possibility  of  an  all  optical  AM-FM 
conversion. 
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Finally,  the  effects  of  group  velocity  splitting  in  coastal  shallow  water 
acoustics  should  be  studied.  For  example,  potential  applications  exist  in  the  use 
of  pulse  splitting  as  a  littoral  environment  mine  countermeasure. 

3.  FUTURE  WORK 

Future  work  should  include  continued  investigation  of  the  pulse  splitting 
phenomena  using  input  signals  of  smooth  modulation.  The  splitting  appears  to 
be  more  easily  observable  for  pulses  of  this  type. 

Future  experiments  should  also  include  determination  of  the  relationship 
between  the  relative  phase  and  the  square  of  the  wave  amplitude  for  a  range  of 
frequencies  above  cutoff.  As  predicted  by  the  theory  presented  in  Chapter  4,  the 
nonlinear  coefficient  changes  sign  for  frequencies  above  1.897  times  the  cutoff 
frequency.  For  this  frequency  range,  modulational  instabilities  lead  to  the 
formation  of  solitons.  To  our  knowledge,  no  experimental  observation  of 
acoustic  solitons  has  been  made.  A  future  experiment  should  explore  this 
possibility. 

Experimental  verification  of  the  predicted  AM-FM  conversion  effect  should 
also  be  pursued. 
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